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An overview is given on how superconductivity with anisotropic pair- 
ing can be realised from repulsive electron-electron interaction, (i) We start 
from the physics in one dimension, where the Tomonaga-Luttinger theory 
predicts that, while there is no superconducting phase for the repulsive case 
for a single chain, the phase does exists in ladders with the number of legs 
equal to or greater than two, as shown both by analytically (renormalisation) 
and numerically (quantum Monte Carlo), (ii) We then show how this pair- 
ing has a natural extension to the two-dimensional case, where anisotropic 
(usually d) pairing superconductivity arises mediated by spin fluctuations 
(usually antiferromagnetic), as shown both by analytically (renormalisation) 
and numerically (quantum Monte Carlo) . (iii) We finally discuss how the su- 
perconductivity from the electron repulsion can be "optimised" (i.e., how Tq 
can be raised) in 2D and 3D, where we propose that the anisotropic pairing 
is much favoured in systems having disconnected Fermi surfaces where Tq 
can be almost an order of magnitude higher. 

1 Introduction 

There is a growing realisation that the high-Tc superconductivity found in the 
cuprates in the 1980's has an electronic mechanism — namely, anisotropic 
pairing from the repulsive electron-electron interaction. Superconductivity 
from electron repulsion is conceptually interesting in its own right, and has 
indeed a long history of discussion. In fact, in the field of electron gas, i.e., 
electron system with the Coulombic electron-electron interaction, Kohn and 
Luttinger 1 pointed out, as early as in the 1960's, that the electron gas should 
become superconducting with anisotropic pairing (having nonzero relative 
angular momenta) at sufficiently low temperatures in a perturbation theory. 
This becomes an exact statement for dilute enough electron gas, where p- 
wave (with the relative angular momentum = 1) should arise, as far as the 
static interaction is concerned |2I3| . 

While these have to do with the long-range Coulomb interaction where 
the dominant fluctuation is charge fluctuation, the problem we would like 
to address here is the opposite limit of short-range repulsion, as appropriate 
for strongly-correlated systems such as transition metal oxides. There, the 
dominant fluctuation is the spin fluctuation. The most widely used model is 
the Hubbard model having the on-site repulsion, U. If the one-band Hubbard 



2 Hideo Aoki 



model, the simplest possible model for repulsively correlated electron systems, 
superconducts, the interest is not only generic but may be practical as well, 
which has indeed been a challenge in the physics of high Tq superconductivity. 

To develop a theory for that, it is instructive to start with one-dimensional 
(ID) systems. When the system is purely ID, we have an exact effective the- 
ory, which is the Tomonaga-Luttinger theory and is exactly solvable in terms 
of the bosonisation and renormalisation. So we start with this, where no 
superconducting phase is shown to exist when the interaction is repulsive. 
When there are more than one chains, which is called ladders, superconduct- 
ing phase appears. If one closely looks at the pairing wavefunction, this is a 
pairing having opposite signs across two bands where the key process is the 
interband pair hopping. 

We then show that this physics has a very natural extension to two- 
dimensional (2D) systems. There, anisotropic (usually d having the relative 
angular momentum of 2) pairing superconductivity can arise. If one looks 
at the pairing wavefunction, this is a pairing having opposite signs across 
the key interband pair-hopping processes. The key process is dictated by the 
peak in the spin structure (usually antiferromagnetic) . 

We finally look at how this kind of anisotropic pairing superconductivity 
can be "optimised" , namely how we can make Tq higher. We first note that 
u Tc is very low in the electron mechanism" in that Tq is usually two orders 
of magnitude lower than the electronic energy. The main reason is the node 
in the gap function, which has to exist for the anisotropic pairing, intersects 
the Fermi surface. So we can propose, and show, that systems that have 
disconnected Fermi surface has much higher Tq- 

2 ID — Tomonaga-Luttinger theory and the physics 
of ladders 

2.1 Tomonaga-Luttinger theory 

It was Tomonaga who pioneered the many-body physics in ID. In his 1950 
paperjl] the essence of the whole idea is already there, although the theory 
is now often called Tomonaga-Luttinger. When the system is ID, the Fermi 
energy, Ep, intersects the band at two points, left-moving branch (L) and 
the right-moving one (R; Fig^a)). The dispersion around these points may 
be approximated as linear functions of the wavenumber, k. When we do this, 
every electron-hole excitation across Ep becomes a creation operator of a 
sound wave (which is a boson) . 

As for the electron-electron interaction, the matrix elements may be clas- 
sified into four categories: (i) backward scattering, where one electron at R 
jumps to L while another from L to R, (ii) forward scattering, where R 
jumps to R, L jumps to L, (iii) umklapp scattering, where {[R,R] jumps 
to [L,L] or vice versa), (iv) forward scattering within each branch ([R, R] 
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Fig. 1. (a) Tomonaga-Luttinger model, in which the low-energy excitations around 
the Fermi energy (shaded) at k = ±&f are considered for ID systems, (b) Weak- 
coupling result for the phase diagram against the on-site repulsion, U, and the 
off-site repulsion, V, for the ID extended Hubbard model at half filling. SS: spin- 
singlet superconductivity, TS: triplet superconductivity. 



to [R,R] or [L. L] to [L, L\). Tomonaga-Luttinger theory |5I7I8I9I10I1T] is 

a weak-coupling theory (i.e., theory for the case when the electron-electron 
interaction is weak enough), where only question low-energy processes. For 
that we can integrate out the higher-energy processes in the perturbational 
renormalisation-group sense. We can then look at the flow of the renormal- 
isation equation, and its end point called the fixed point. To discuss the 
nature of the fixed-point Hamiltonian, it is convenient to bosonise (i.e., to 
write everything in terms of boson operators). The final result for the ef- 
fective Hamiltonian is written in terms of two boson fields, spin phase (<j)) 
and charge phase (9), whose stiffness (coefficients of (dtp) 2 , (d9) 2 ) is given in 
terms of only two quantities, K a , K p , which determine everything, including 
whether the ground state is superconducting. To be more precise, in ID even 
a "long-range" order can only have a two-point correlation that decays with 
a power law (cx l/r") where the exponent a is dictated, for each of the order 
parameters considered, by K ai K p . 

For every Hamiltonian originally given, we can calculate the four scatter- 
ing parameters, and then renormalise them. If we look at the phase diagram 
(Fig^b)) for the extended Hubbard model (where we have an off-site inter- 
action, V, on top of the on-site one, U), there is no superconducting phase 
when all the interactions U and V are repulsive (> 0). 
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Incidentally, there is no magnetism, either, for a single chain. This is due 
to the well-known Lieb-Mattis theorem, which dictates that electrons in ID 
are entirely non-ferromagnetic. The proof makes use of the fermion statistics 
of electrons, where a key factor is no two electrons can pass each other in ID, 
or, in the words in Mattis's textbook^, "neighbours remain neighbours till 
death did them part" . 

2.2 Pairing in ladders, or 1 + 1 ^ 2 

When there are more than one chain with inter-chain hopping and/or inter- 
action, the physics can be, and is indeed, entirely different. The model, then, 
becomes multi-band (i.e., n-band system for n-leg ladder). The Fermi energy 
can intersect the dispersion at 2n points (Fig|2J|, so the model is what can 
be called multiband Tomonaga-Luttinger model. The multiband Tomonaga- 
Luttinger model has been studied in various context, including the excitonic 
phase in electron- hole systems |12|. transport properties and interband ex- 
citations in quantum wires as detected by Raman spectroscopy by Sassetti 
et al.[H| 




Fig. 2. Multi-band ID model (right panel) for the Hubbard model on a ladder 
(left). An oval on the ladder represents the inter-chain pairing. 
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In the context of the high Tq, the idea of superconductivity in multi-chain 
(or "ladder") systems was kicked off theoretically in 1986, when Schulz 15 
proposed a possible relation between ladders at half-filling and Haldane's 
conjecture for spin chains. He made the following reasoning: If we consider 
repulsively interacting electrons on a ladder, the undoped system will be a 
Mott insulator, so that we may consider the system as an S = 1/2 antiferro- 
magnetic (AF) Heisenberg magnet on a ladder for large Hubbard repulsion U. 
Schulz 's analysis is that an AF S — n/2 single chain, which is exactly the 
Haldane's system |16I17| . is similar to an S — 1/2 AF ladder with n-legs. For 
the spin chains, Haldane ^B] has conjectured that the spin excitation should 
be gapless for half-odd-integer spins {n: odd) or gapful for integer spins (n: 
even). If the situation is similar in ladders, a ladder having an even num- 
ber of legs will have a spin gap, associated with a 'spin-liquid' ground state 
where the quantum fluctuation is so large that the AF correlation decays 
exponentially 

Dagotto et aL^H] and by Rice et a/.^Hji then suggested the possibility of 
superconductivity associated with the spin gap. The presence of a spin gap, 
i.e., a gap in the spin excitation which is indicative of a quantum spin liquid, 
in the two-leg ladder (or, more generally, in even legs) is a good news for 
superconductivity, since an idea proposed by Anderson |20| in the context of 
the high-Tc superconductivity suggests that a way to obtain superconductiv- 
ity is to carrier-dope spin-gapped systems. The superconductivity in even-leg 
ladders is in accord with this. Subsequently superconductivity has been re- 
ported for a cuprate with a ladder structure |21j. although it later turned out 
that this material has a rather strong two-dimensionality that may dominate 
the superconductivity. 

For doped systems, the conjecture for superconductivity 19 is partly 
based on an exact diagonalisation study for finite t-J model on a two-leg lad- 
der. [3U| This was then followed by analytical and numerical 32 33 34I35| 
works on the doped t-J ladder, for which the region for the dominant pairing 
correlation appears at lower side of the exchange coupling J than in the case 
of a single chain. 

On the other hand the Hubbard model on a ladder is of general interest. |3*rJ| 
Although the Hubbard crosses over to the t-J model for U — > oo, we have 
only an infinitesimal J there, so the result for t-J model does not directly 
answers this. Since there is no exact solution for the Hubbard ladder, we can 
proceed in two ways: for small U we can adopt an analytic method, which 
is the weak-coupling renormalisation-group theory, where the band structure 
around the Fermi points is linearised in the continuum limit to treat the inter- 
action with a perturbative renormalisation group. The weak-coupling theory 
with the bosonisation and renormalisation-group techniques has been applied 
to the two-leg Hubbard ladder. j22l23l24l25<2"6T27] 
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The Hamiltonian of the two-leg Hubbard ladder is given in standard no- 
tations as 

« = -* E + h - c -) - ^ E( c - c - + h - c -) 

a{i)u icr 

+uY f nf t nf l , (1) 

ai 

where a(= 1, 2) specifies the chains, or in the momentum space as 
H = -2t 5^ - 2t ± £ 4UI 

fika k(7 

+U ^^(interaction of the form c'c^cc), (2) 

where \i specifies the bonding (/i = 0) and anti-bonding (/i = tt) bands, so 
labelled since k y = 0, n, respectively. 

The part of the Hamiltonian, Hd, that can be diagonalised in the bosoni- 
sation includes only intra- and inter-band forward-scattering processes arising 
from the intrachain forward-scattering terms. We can then define bosonic op- 
erators as in the single-chain case. If we introduce the phase variables as in 
the single-chain case, Hd, written in terms of them, is separated into the spin- 
part -ffspm and the charge-part -ff c h a rgc- While H sp i n is already diagonalised, 
^charge can be made so with a linear transformation, and the diagonalised 
H charge is written in terms of the correlation exponent K p i,i = 1,2. So we 
end up with the total Hamiltonian that reads 

H = kineticenergy + Hd + pair-hopping terms. (3) 

Here, the pair-hopping (or pair-scattering) term represents those part of the 
interaction Hamiltonian, in which a pair of electrons is scattered via the 
interaction to another pair of electrons. 

At half-filling, the system reduces to a spin-liquid insulator having both 
charge and spin gaps When carriers are doped to the two-leg Hubbard 
ladder, on the other hand, the relevant scattering processes at the fixed point 
in the renormalisation-group flow are the pair hopping across the bonding and 
anti-bonding bands (FigEJ, c^c^c^c® + h.c. in k space, and the backward- 
scattering process within each band. The importance of the pair-hopping 
across the two bands for the dominance of pairing correlation in the two- 
leg Hubbard ladder is reminiscent of the Suhl-Kondo mechanism, which was 
proposed back in the 1950's for superconductivity in a quite different context 
of the s-d model for the transition metals. |37I38| 

The renormalisation results in a formation of gaps in both of the two spin 
modes and a gap in one of the charge modes. This leaves one charge mode 
massless, where the mode is characterised by a critical exponent K p . Then 
the correlation of the intraband singlet pairing, 

WcLeV* - c w c -fc,- J, (4) 
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decays like 1 / r 1 / ( 2K p\ where K p should be close to unity in the weak-coupling 
regime. So this should be the dominant phase, which is, expressed in real space 
as c\ a cf _ a — c] _ a .cf rT , an interchain singlet pairing. 

2.3 How to detect pairing in quantum Monte Carlo studies? 

The perturbational renormalisation group is in principle guaranteed to be 
valid only for sufficiently small interaction strengths (U <C t), so that its 
validity for finite U (~ t) has to be checked. To be more precise, the renor- 
malisation approach can tell whether the interaction flows into weak coupling 
(with the relevant mode gapless) or into strong-coupling regime (gapful) for 
small enough interactions, but the framework itself (i.e., the perturbational 
expansion) might fail for stronger interactions. 

This is where numerical studies come in. Numerical calculations for finite 
U have been performed with the exact diagonalisation, DMRG or quantum 
Monte Carlo (QMC) methods, |39I4()I41I42I43I44| but in an earlier stage the 
results are scattered, where some of the results seemed inconsistent with 
the weak-coupling prediction: a DMRG study by Noack et al. for the doped 
Hubbard ladder shows the enhancement of the pairing correlation over the 
U = result strongly depends on the inter-chain hopping, t± 28 45 . Quan- 
tum Monte Carlo (QMC) results also exhibit an absence|3Bl or presence |4*7| of 
the enhancement depending on the hopping parameters and/or band filling. 

Recently, however, a QMC study by Kuroki et al has resolved the 
puzzle, and has clearly detected an enhanced pairing correlation. A key fac- 
tor found there in detecting superconductivity in any numerical calculation, 
which also resolves the origin of the former discrepancies is: we have to ques- 
tion a very tiny energy scale (<C starting electronic energy scale, t, U) in 
detecting the pairing. This immediately implies that the discreteness of en- 
ergy levels in finite systems examined in QMC studies enormously affects 
the pairing correlation — If the level separation is greater than the energy 
scale we want to look at, any feature in the correlation function will be easily 
washed out. This can be circumvented if we tune the parameters so as to 
make the separation between the levels just below and above Ep tiny (i.e., 
to make the LUMO-HOMO nearly degenerate in the quantum chemical lan- 
guage), which should be a reasonable way to approach the bulk limit where 
the levels are dense. The importance of small offsets between the highest oc- 
cupied and lowest unoccupied levels has also been stressed by Yamaji et al. 
for small systems. 03 

So we have applied the (projector) Monte Carlo method^Hl to look into 
the ground state correlation function P(r) = (Oj +r Oi) of this pairing for 
finite values of Z7(~ t). We show in Fig|2|the result for P(r) for tj_ = 0.98, 
U = 1 and the band filling n = 0.867 = 52 electrons/ (30 rungs x 2 sites). 
The U = result (dashed line) for these two values of t± are identical because 
the Fermi sea remains unchanged. If we turn on U, we can see that a large 
enhancement over the U = result emerges at large distances. 
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Fig. 3. The pairing correlation function plotted against the real space distance 
r in a 30-rung Hubbard ladder having 56 electrons for U — 1 with t± = 1.975 
(square) .|48| The dashed line is the non-interacting result. 

We have deliberately chosen the value of t± = 0.98 to make the one- 
electron energy levels of the bonding and anti-bonding bands lie close to 
each other around the Fermi level within 0.004i. This is much smaller than 
the energy scale we question (which is the spin gap in the present ID case). 
In fact, a 5% change in t ± = 0.98 -> 1.03, for which the LUMO-HOMO 
separation blows up to ~ O.lt, washes out the enhancement in the correlation 
function. In the latter case the renormalisation of higher energy modes has 
to stop at this energy scale, so that the interband pair hopping process will 
not be renormaliscd into a strong coupling, while in the weak-coupling theory 
the renormalisation all the way down to the Fermi level is assumed. 

While it is difficult to determine the decay exponent of the pairing correla- 
tion .P(r), we can fit the data by assuming a trial function expected from the 
weak-coupling theory, P(r) oc c/\fr + (1 — c)/r 2 + [cos(2fc^r) + cos(2fc J.r)]/r 2 
where the overall decay at large distances is assumed to be oc 1/ y/r as dictated 
in the weak-coupling theory. This form reproduces the result surprisingly ac- 
curately. 
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3 Three-leg ladder and 1D-2D crossover 

Now, the physics of ladders can provide quite an instructive line of approach 
for understanding the physics in two-dimensional systems via the crossover 
from ID to 2D (two-dimensions). So let us first look at the three- leg ladder. 

3.1 Three-leg ladder 

If we return to ladders, one can naively expect that ladders with odd-number 
(e.g., 3) of legs will have no spin gap, which would then signify an absence of 
dominating pairing correlation ('even-odd conjecture for superconductivity'). 
As far as the spin gap in undoped ladders is concerned, experiments on a 
class of cuprates, Sr„_iCu n 02 n -i having n-leg ladders, have supported the 
conjecture. |51I52I53I54<55] So it was believed that odd-numbered legs only 
have the usual 2kp spin-density wave (SDW) rather than superconductivity. 

Kimura et al |56!57j .however, showed that that is too simplistic a view, and 
that, while the even-odd conjecture for the spin gap is certainly correct, an 
odd-number of legs does indeed superconduct by exploiting the spin-gapped 
mode. In that work the pairing correlation in the three-leg Hubbard ladder 
has been examined |5()|. 

We start with the weak-coupling theory for correlation functions for the 
three-leg Hubbard ladder (Fig0J . Arrigoni has looked into a three- leg ladder 
with weak Hubbard-type interactions with the perturbational renormalisation- 
group technique to conclude that gapless and gapful spin excitations coexist 
in three legs. [HE] 

He has actually enumerated the numbers of gapless charge and spin modes 
on the phase diagram spanned by the doping level and the interchain hopping, 
t±. He found that, at half-filling, one gapless spin mode exists. For general 
band filling, one gapless spin mode remains in the region where the Fermi 
level intersects all the three bands in the noninteracting case. From this, 
Arrigoni argues that the 2k p SDW correlation should decay as a power law 
as expected from experiments. Arrigoni's result indicates that two gapful spin 
modes exist in addition. The charge modes, on the other hand, consists of 
two gapless modes and one gapful mode. 

The question we address then is what happens when gapless and gapful 
spin modes coexist. This is an intriguing problem, since it may well be pos- 
sible that the presence of gap(s) in some out of multiple spin modes may be 
sufficient for the dominance of a pairing correlation. Schulz |59j has indepen- 
dently shown similar results for a subdominant 2k p SDW and the interchain 
pairing correlations. 

Physically, the picture that emerges as we shall describe below, is that 
the two spin gaps, which are relevant to the pairing, arise as an effect of 
the pair-hopping process that is the many-body matrix element transferring 
two electrons simultaneously across the outermost bands (i.e., the top and 
bottom bands for a three- leg ladder) (Fig^J. In this sense the mechanism is 
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Fig. 4. Three-leg ladder (inset) and its dispersion. The dashed arrows represent 
the pair-hopping process, and ovals represent the inter-chain pairing. 



reminiscent of the situation in the two-leg case or the Suhl-Kondo mecha- 
nism. |37I38| 

The correlation functions can be calculated with the bosonisation method [H] 
for the three-leg Hubbard model. 

We can define three bosonic operators, where we diagonalise -ff c hargo in 
terms of the three correlation exponents K pi , i = 1,2,3. As Arrigoni pointed 
out [SHj the pair-hopping processes across the top and bottom bands become 
relevant as the renormalisation is performed. In order to actually calculate the 
correlation functions, we have to express the relevant scattering processes in 
terms of the phase variables. The fixed-point Hamiltonian density, H* , takes 
the form, in terms of the phase variables, 

H* oc - 3back(l)cos[2(/) 1+ (2;)] - g bac k(3)cos[2(/> 3+ (iz;)] 

+ 2gp airho pp ing (l,3)cos[V2x 1 ~(x)]sin(/) 1+ (j:)sin(/i 3+ (x), (5) 

where gback(l), <7back(3) are negative large quantities, and g p h(l,3) is a posi- 
tive large quantity. This indicates the following. Two spin phases, </> 1+ , </> 3+ , 
become long-range ordered and fixed, respectively, while 4> 2+ is not fixed to 
give a gapless spin mode. Similarly, the difference in the charge phases, 9 l , 
for the outermost bands, 



1 



(6) 



Superconductivity from the repulsive electron interaction 11 



is ordered and fixed, and the charge gap opens for this particular mode. 

Now we can calculate the correlation functions, since the gapless fields 
have already been diagonalised, while the remaining gapful fields have the 
respective expectation values. 

Among various order parameters, the dominant one (with the longest tail 
in the correlation) is the singlet pairing across the central and edge chains 
(Fig0J, which is, in the band picture, 

o d ~ - (7) 

a 

We call this "d" for the following reason. Since we have taken a continuum 
limit along the chain, it is not straightforward to name the symmetry of a 
pairing. However, we could call the above pair as d-wave-like in that the 
pairing, in addition to being off-site on the rung, is a linear combination of 
a bonding band and an anti-bonding band with opposite signs. Since the 
relevant pair-hopping is across these bands, we can say that there is a node 
in the pair wavefunction along the line that bisects the relevant pair-hopping. 
Calculation of the correlation function of Od gives 

(O d (x^O d (0)) ~ x ~ Hl ^ + ^\ (8) 

which is the dominant ordering. In the weak-interaction limit (U — > +0), 
where all the i*sT*'s tend to unity, the "d" -pairing correlation decays as slowly 
as x _1 / 2 , while the other correlations decay like x~ 2 . We can see that the 
interchain pairing exploits the charge gap and the spin gaps to reduce the 
exponent of the correlation function, in contrast to the intrachain pairing. 
Namely, we would have (—1) in the exponent if the spin were gapless. This 
alone would only result in a 1/r decay, but the charge mode xi- is further 
locked, which further reduces the exponent (down to \j\fr in the limit U — ► 

Now, how the pairing correlation in the three-leg Hubbard ladder looks 
like when U is finite? Our QMC result for the three-leg Hubbard ladder 
exhibits an enhancement of the pairing correlation even for finite coupling 
constants, U/t = 1 ~ 2. j^ZI As in the two-leg case with a finite U, we have 
taken care that levels below and above the Fermi level are close. 

3.2 1D-2D crossover 

We have seen that the weak-coupling theory (pcrturbational renormalisation 
+ bosonisation) predicts that the interband pair hopping between the in- 
nermost and outermost Fermi points (k x ,k y ) « (±fc^, 0), (±fc£, tt) becomes 
relevant (i.e., increases with the renormalisation). This concomitantly makes 
the two-point correlation of the interchain singlet decay slowly with distance. 
In fc-space the dominant component of this pair reads 
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Now, when Ep intersects the outermost-band top and the innermost-band 
bottom with k F ~ tt, k F ~ (Fig[5]left), intrachain nearest-neighbor singlet 
pair also has a dominant Fourier-component equal to eqn.JHJ with a phase 
shift 7r relative to the interchain pairing. Thus, a linear combination which 
amounts to the d x 2_ y 2 pairing should become dominant. We shall sec this is 
exactly what happens in the 2D squre lattice around the half filling, Fig[5] 
right. 




Fig. 5. Band dispersion of the square lattice, on which the Fermi surface for Ef 
close to half filling and a typical pair-hopping process (dashed arrows) are shown. 
Similar plot for a ladder is attached for comparison. 



4 Superconductivity from the repulsive interaction in 
2D 

As we have seen, the two-leg and three-leg Hubbard ladders do superconduct, 
so what will happen if we consider n-leg ladders for n = 3,4, ...,oo to reach 
the 2D square lattice. This view enables us to have a fresh look at the 2D 
Hubbard model. 

So we move on to the repulsive Hubbard model on the square lattice. 
The seminal notion that the high- Tq conductivity in cuprates should be re- 
lated the strong electron correlation was first put forward by Anderson. |61| 
There the superconductivity is expected to arise from the pairing interaction 
mediated by spin fluctuations (usually antiferromagnetic). A phenomenology 
along this line such as the self-consistent renormalisation 62 63 64 65 has 
succeeded in reproducing anisotropic <i-wave superconductivity. Microscopi- 
cally, the repulsive Hubbard model, a simplest possible model for correlated 
electrons, should capture the physics in cuprates. |60) Some analytical calcu- 
lations have suggested the occurrence of d ;E 2_ y 2-wave superconductivity in 
the 2D Hubbard model. |66I67I68I69I7U| In particular, fluctuation exchange 
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approximation (FLEX) , developed by Bickers et at [7T| , has also been applied 
to the Hubbard model on the square lattice |72I73| to show the occurrence 
of the superconductivity. 

Numerical calculations have also been performed extensively. Finite bind- 
ing energy [74175) and pairing interaction vertex [?6l77l78| were found in those 
calculations. Variational Monte Carlo calculations show that a superconduct- 
ing order lowers the variational energy. 79 80 Nevertheless, there had been 
a reservation against the occurrence of superconductivity in the Hubbard 
model because the pairing correlation functions do not show any symptom 
of long-range behavior in some of the works. [?8I81I82| 

Again, Kuroki et al0H] showed for the first time that QMC does indeed 
exhibit symptoms of superconductivity if we take proper care of a small 
energy scale involved, i.e., the d-wave pairing correlation becomes long-tailed 
when the Fermi level lies between a narrowly separated levels residing on the 
k-points across which the dominant pair hopping occurs. An enhancement 
of the pairing correlation has in fact been found by exact diagonalisationJJJ] 
and by density matrix renormalisation group [43] when Ep lies close to the 
fc-points (0, it) and (it, 0). Although the d a ,2_ l) 2-like nature of the pairing was 
suggested, jUjj d x 2_ y 2 pairing correlation itself has not been calculated. So, in 
our quest for 2D, we first calculate the correlation function with QMC. Here 
we employ the ground-state, canonical-ensemble QMC,jHH| where we take the 
free Fermi sea as the trial state. 



4.1 Anisotropic pairing in 2D 

First, let us look at why the attractive interaction is by no means a nec- 
essary condition for superconductivity, which can quite generally arise from 
repulsive electron-electron interactions, which seems to be still not realised 
well enough. If we look at the BCS gap equation, we can immediately see 
that superconductivity can readily arise from repulsive interactions. The gap 
equation reads 



1 = -^E7777T tanh \lfe( k ') 



k 2e(fe,) 



2" 



(V(k,k')A(k)A{k')) FS 

v * = — wm^s — ' (10) 

where s(k) is the band energy measured from the chemical potential, f3 — 
l/(fcsTc), V(k,k) the pair-hopping matrix element, A(k) the BCS gap 
function, and (...)fs is the average over the Fermi surface. So, if A(k) has 
nodes across k <-> k' (i.e., changes sign before and after the pair hopping), 
the originally repulsive V > acts effectively as an attraction, 



(V(k,k')A(k)A(k')) FS < 0. 
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This most typically happens for the d x 2_ y 2 pairing with A(k) oc cos(k x ) — 
cos(k y ) when the dominant pair hopping occurs across k ~ (0,tt) «-> fe' ~ 

(7T,0). 

When the spin fluctuation is antiferromagnetic, most typically in bipartite 
lattices such as a square lattice, the importance of the interactions around 
(0, 7r) and (tt, 0) in the 2D Hubbard model has been suggested by various 
authors. |(iVlH8l(i9lV()IVVl8()l8,^l8(il8Vl88| 

Group theoretically, the square lattice has a tetragonal symmetry, so that 
everything, including the gap function, should be an irreducible represen- 
tation of the tetragonal group. The d x 2_ y 2 pairing indeed belongs to Bi g 
representation of this group. 

4.2 Quantum Monte Carlo study for the 2D Hubbard model 

In the context of our QMC study for the 2D Hubbard model, we have to 
take finite systems that have the fc-points around (0, tt) and (tt, 0) close in 
energy. Nameley, our expectation from the study on ladders is that the pair 
hopping processes across around (0, tt) and (tt, 0) may result in d x 2_ y 2 pairing, 
^^[cos(fc x ) — cos(k y )]cf : ^c_f ( .^ in 2D, but an enhanced pairing correlation 
should be detected only when the level offset between the discrete levels 
around those points is small. 

We take 78 electrons in 10 x 10 sites (n = 0.78) with t y = 0.999 with 
periodic boundary condition in both directions. We have taken t y = 0.999, 
because the number of electrons considered here would have an open shell 
(with a degeneracy in the free-electron Fermi sea) for t y = 1, which will 
destabilise QMC convergence. Taking t y — 0.999 lifts the degeneracy to give 
a tiny (< 0.01) but finite Ae°. In FigElwe plot the d x 2_ y 2 pairing correlation, 
defined as 

P(r) = Yl (0\x + Ax,y + Ay)0(x,y)), 

\Ax\ + \Ay\=r 

^(•^5 V) — ^ &(Cx,y,aCx+8,y, — <T ^x,y,(7^x,y-\-S,—a)i (H) 

«=±1,<T 

where the correlation for U — 1 is clearly seen to be enhanced over that for 
U = especially at large distances. 

We can readily show that when the level spacing becomes too large (e.g., 
Ae° ~ 0.1), the enhancement is washed out. In the present choice the energy 
levels around (0,tt), (tt,0) are close (< O.Olt), while the other levels lie more 
than ~ O.lt away from Ep. One might thus raise a criticism that the scat- 
tering processes involving the states away from E F are unduly neglected. We 
can however show (not displayed here) that when other levels exist around 
Ep an enhanced d x 2_ y 2 correlation is obtained as well. 

How about the band-filling (n) dependence? We have calculated the long- 
range part of the correlation, S = J2 r >3 ^ > ( r )' f° r various values of n keeping 
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Fig. 6. QMC result for the d x 2_ y 2 pairing correlation for a 10 x 10 square lattice 
with 78 electrons for U = 1 and t y = 0.999 (square ~).|48| The dashed line represents 
the noninteracting case. Inset depicts the square lattice with the ovals representing 
a d-wave pair. 

Ae a < O.Olt throughout. The result, displayed in Fig0 shows that the en- 
hancement in S for U — 1 has a maximum around a finite doping. Thus the 
message here is that the d x 2_ y 2 pairing is favoured near, but not exactly at, 
half-filling. The fact that the better nesting does not necesarrily imply the 
more enhanced pairing correlation has also been shown in another numerical 
work in the contex of an organic superconductivity. 89 

5 Which is more favourable for superconductivity, 2D 
or 3D? 

The theoretical results described so far indicate that the superconductivity 
near the AF instability in 2D has a 'low Tc ~ O(0. Oli) (t: transfer integral), 
i.e., two orders of magnitude smaller than the original electronic energy, but 
still 'high T c ' ~ O(100 K) for t ~ 0(1 eV). So, identifying the conditions 
for higher Tc in the repulsive Hubbard model is one of the most fascinating 
goals of theoretical studies. For instance, while the high-Tc cuprates are 
layer-type materials with Cu2 planes in which the supercurrent flows, the 
question is whether the two-dimensionality is promoting or degrading the 
superconductivity. 

Hence Arita et al[H0] have questioned: (i) Is 2D system more favourable for 
spin-fluctuation mediated superconductivity than in three dimensions(3D)? 
(ii) Can other pairing, such as a triplet p-pairing in the presence of ferro- 
magnetic spin fluctuations, become competitive? We take the single-band, 



16 Hideo Aoki 



CO 




Fig. 7. The integrated pairing correlation plotted against the band filling n for a 
10 x 10 lattice with U = (o) or U = 1 (square). EH] 



repulsive Hubbard model as a simplest possible model, and look into the 
pairing with the FLEX method in ordinary (i.e., square, trianglar, fee, bee, 
etc) lattices in 2D and 3D. The FLEX method has an advantage that systems 
having large spin fluctuations can be handled. 

As for 3D systems, Scalapino et aZ[5J showed for the Hubbard model 
that paramagnon exchange near a spin-density wave instability gives rise to 
a strong singlet <i-wave pairing interaction, but Tq was not discussed there. 
Nakamura et al 92 extended Moriya's spin fluctuation theory of superconductivity 63 
to 3D systems, and concluded that Tc is similar between the 2D and 3D 
cases provided that common parameter values (scaled by the band width) 
are taken. However, the parameters there are phenomelogical ones, so we 
wish to see whether the result remains valid for microscopic models. 

As for the triplet pairing, the possibility of triplet pairing mediated by fer- 
romagnetic fluctuations has been investigated for superfluid 3 He|2S], a heavy 
fermion system UPtal^J, and most recently, an oxide Sr 2 Ru04[HSl- It was 
shown that ferromagnetic fluctuations favour triplet pairing first by Layzer 
and Fay|2] before the experimental observation of p-wave pairing in 3 He. 
For the electron gas model, Fay and Layzerj^j and later Chubukov 96] has 
extended the Kohn-Luttinger theorem^ to p-pairing for 2D and 3D elec- 
tron gas in the dilute limit. Takada[9Tj discussed the possibility of p-wave 
superconductivity in the dilute electron gas with the Kukkonen-Overhauser 
model 98 . As for lattice systems, 2D Hubbard model with large enough next- 
nearest-neighbor hopping (t') has been shown to exhibit p-pairing for small 
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band fillings. |M| Hlubina |lU(J| reached a similar conclusion by evaluating the 
superconducting vertex in a perturbative way. |101| However, the energy scale 
of the p-pairing in the Hubbard model, i.e., Tc, has not been evaluated so 
far. 

Here we show that (i) d-wave instability mediated by AF spin fluctuation 
in 2D square lattice is much stronger than those in 3D, while (ii) p-wave in- 
stability mediated by ferromagnetic spin fluctuations in 2D are much weaker 
than the d-instability. These results, which cannot be predicted a priori, sug- 
gest that for the Hubbard model the 'best' situation for the pairing instability 
is the 2D case with dominant AF fluctuations. 

We consider the single-band Hubbard model with the transfer energy 
tij = t(= 1 hereafter) for nearest neighbors along with tij — t' for second- 
nearest neighbors, which is included to incorporate the band structure depen- 
dence. The FLEX starts from a set of skeleton diagrams for the Luttinger- 
Ward functional to generate a (fc-dependent) self energy based on the idea of 
Baym and Kadanoff 102,. Hence the FLEX approximation is a self-consistent 
perturbation approximation with respect to on-site interaction U. 

To obtain Tq, we solve the eigenvalue (Eliashberg) equation, 



XU^(k) = ^E^(k')\G(k')\ 2 V^(k-k'), 



N 



(12) 



where (k) is the anomalous self energy, k = (k, iuj n ) with ui n = (2n—l)nT 
being Matsubara frequencies, and the pairing interaction, V^ 2 \ comprises 
contributions from the transverse spin fluctuations, longitudinal spin fluctu- 
ations and charge fluctuations, namely, 



V {2 \k,k') = -U 2 
1 



:Xch(k - k') 



x zz (k-k') + x ± {k + k') 



1 - U 2 xUk - k>) 



U 2 XiAk + k') 

1 - U X irr(k + k>) 



Here Xch is the charge susceptibility, X zz (x ± ) the longitudinal (transverse) 
spin susceptibility, and 



X 



irr (g) = -(T/7V)^G(fc)G(fc + g ) 



the irreducible susceptibility constructed from the dressed Green's function. 
The dressed Green's function, G(k), obeys the Dyson equation, 



Giky 1 =G°(k) 1 -S(k), 



(13) 
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where G° is the bare Green's function, and £ the self energy with 

^( fc ) = ^E G ( fc -^ (1) (<z)- (14) 



If we take RPA-type bubble and ladder diagrams for the interaction 
have 



wo 



V {1) (q) = \u 2 X irM 



1 



l + C^Xirr(g) 
1 



1 - UXirriq) 



which completes the set of equations. 

Since we have £^ 2 \k) = £( 2 \~k) for the spin-singlet pairing whereas 
£( 2 \k) = —£( 2 '(—k) for the spin-triplet pairing, V^ 2 \k, k') becomes a func- 
tion of k — k' = q with 



V^{q) = -\ 
for the singlet pairing, and 

V^{q) = \ 



U 2 X irr(q) 



i - UxUq) 



U 2 XirM 



U 2 X irr(q) 



1 - UXirr(q)_ 



i + UxivM 



^Xirr(g) 

i + UxiAq) 



(15) 



(16) 



for the triplet pairing. T — Tc is identified as the temperature at which the 
maximum eigenvalue AMax reaches unity. 

Let us start with the 2D case having strong AF fluctuations. We have 
first obtained xrpa(<?) = Xo/0- ~ Uxo) as a function of the momentum for 
the Hubbard model on a nearly half-filled (n = 0.85) square lattice, where a 
dominant AF spin fluctuation is seen as Xrpa peaked around (jr, it). We can 
then plug this into the Eliashberg equation to plot in Fig[H{a) Aj^ax as 
a function of temperature T. 

Tc is identified as the temperature at which AMax becomes unity, which 
occurs at T ~ 0.02 for the square lattice, in accord with previous results [721103] . 

If we move on to the case with ferromagnetic spin fluctuations where 
triplet pairing is expected, this situation can be realised for relatively large 
t'(~ 0.5) and away from half-filling in the 2D Hubbard model. Physically, the 
van Hove singularity shifts toward the band bottom with t' , and the large 
density of states at the Fermi level for the dilute case favours the ferromag- 
netism. We have found that AMax becomes largest for n — 0.3, t' = 0.5. 

Xrpa is indeed peaked at r (k = (0, 0)). The question then is the behavior 
of AMax as a function of T, which shows that AMax is much smaller than that 
in the AF case. 

A low Tc for the ferromagnetic case contrasts with a naive expectation 
from the BCS picture, in which the Fermi level located around a peak in 
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Fig. 8. (a) The maximum eigenvalue of the Eliashberg equation against tempera- 
ture for the Hubbard model on a square lattice with n — 0.85 and U = 4 and on 
a cubic lattice with n — 0.8, the second-neighbour hopping t' = —0.2 and U = 8. 
x(k,0) (b) as a function of wavenumber and Imx(feMaxjW) (c) (normalised by its 
maximum) as a function of ui/t at T — 0.03t are also shown for the two lattices. 
Top right panel schematically shows the regions, in respective dimensions, that 
contribute to the pairing. 

the density of states favours superconductivity. We may trace back two-fold 
reasons why this does not apply. First, if we look at the dominant (oc 1/[1 — 
Uxo(q)]) term of the pairing potential V 1 - 2 ^ itself in eqs. I|15|l and l|16() . the 
triplet pairing interaction is only one-third of that for singlet pairing. Second, 
the factor \G\ 2 for the ferromagnetic case is smaller than that in the AF 
case, which implies that the self-energy correction is larger in the former. 
Larger self-energy (smaller \G\ 2 ) works unfavourably for superconductivity 
as seen in the Eliashberg equation (|12|l . When we take a larger repulsion 
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U to increase the triplet pairing attraction (susceptibility), this makes the 
self-energy correction even stronger. 

Let us now move on to the case of d-wave pairing in the 3D Hubbard 
model, for which FLEX was first applied by Arita et allHO]- In simple-cubic 
systems, we find that the r% representation of Oh group[HI] has the largest 
-Wax- We have found that AMax for this symmetry becomes largest for n = 
0.8, t' — —0.2 ~ —0.3 and U — 8 ~ 10. In Fig. IHfa), we superpose AMax as 
a function of T, where we can immediately see that the pairing tendency in 
3D is much weaker than that in 2D. 

Why is the d-superconductivity much stronger in 2D than in 3D? We can 
pinpoint the origin by looking at the various factors involved in the Eliashberg 
equation. Namely we question the height of and \G\ 2 along with the 
width of the region, both in the momentum sector and in the frequency 
sector, over which V^ 2 \k) contributes to the summation over k = (k,iu> n ). 
We found that the maximum of |G| 2 is in fact larger in 3D than in 2D. 
The width of the peak in xrpa on the frequency and momentum axes is 
surprisingly similar between 2D and 3D as displayed in FigJSfb)(c). Note 
that if the frequency spread of the susceptibility scaled not with t but with 
the band width, as Nakamura et aZ 92] have assumed, AMax would have become 
larger. Now, A in the Eliashberg equation 1)12(1 is cx (a/L) D , where L is the 
linear dimension of the system and a the width in the momentum space for 
the effective attraction, this factor is much smaller in 3D than in 2D as far as 
the main contribution of to the pairing occurs through special points in 
the k-space (e.g., (tt, tt) or (n, tt, tt) for the antiferromagnetic spin fluctuation 
exchange pairing). So we can conclude that this is the main reason why 2D 
is more favourable than 3D. 

To summarise this section, d-pairing in 2D is the best situation for the 
repulsion originated (i.e., spin fluctuation mediated) superconductivity in the 
Hubbard model. Monthoux and Lonzarich |104j have also concluded for 2D 
systems, by making use of a phenomenological approach, that the d-w&ve 
pairing is much stronger than p-wave pairing, which is consistent with the 
present result. In this sense, the layer- type cuprates do seem to hit upon the 
right situation. 

This is as far as one-band model having simple Fermi surfaces are con- 
cerned. Indeed, if we turn to heavy fermion superconductors, for instance, in 
which the pairing is thought to be meditated by spin fluctuations, the Tc, 
when normalized by the band width W, is known to be of the order of 0.001PF. 
Since the present result indicates that Tc, normalized by W, is ~ 0.0001VF 
at best in the 3D Hubbard model, we may envisage that the heavy fermion 
system must exploit other factors such as the multiband. Neverthless, recent 
experimental finding [l05j that a heavy-fermion compound Ce(Rh,Ir,Co)In5 
has the higher T c for the more two-dimensional lattice (with larger c/a) is 
consistent with our prediction. 
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6 How to realise higher Tc in anisotropic pairing — 
disconnected Fermi surfaces 

Ironically, the main question about the superconductivity from the electronic 
mechanism is "why is Tc so low?" , which has been repeatedly raised in lit- 
eratures. Namely, one remarkable point is T c ~ O(0. Oli), esitmated for the 
repulsive Hubbard model in the two-dimensional (2D) square lattice, is two 
orders of magnitudes smaller than the starting electronic energy (i.e., the 
hopping integral t), although this gives the right order for the curates' T c . 
We have seen that even the best case, as far as these ordinary lattices are 
concerned, has T c ~ 0(0. Oli). As discussed in Ref. |106| . there are good rea- 
sons why T c is so low: One reason is the effective attraction mediated by spin 
fluctuations is much weaker than the original electron-electron interaction, 
U. Another important reason is the presence of nodes in the superconducting 
gap function greatly reduces T c : While the main pair-scattering, across which 
the gap function has opposite signs to make the effective interaction attrac- 
tive, there are other pair scatterings around the nodes that have negative 
contributions to the effective attraction by connecting fc-points on which the 
gap has the same sign. 

So a next important avenue to explore is: can we improve the situation 
by going over to multiband systems. Kuroki and Arita |l()f>] have shown that 
this is indeed the case if we have disconnected Fermi surfaces. In this case T c 
is dramatically enhanced, because the sign change in the gap function can 
avoid the Fermi pockets, where all the pair-scattering processes contribute 
positively. |10fill07| This has been numerically shown to be the case for the 
triangular lattice (for spin-triplet pairing) 108 and a squre lattice with a 
period-doubling |106| . where T c as estimated with FLEX is as high as O(0. It). 

To be more precise, the key ingredients are: (a) when the Fermi surface is 
nested, the spin susceptibility x(q, u>) has a peak, (b) When a multiband sys- 
tem with a disconnected Fermi surface has an inter-pocket nesting (i.e., strong 
inter-pocket pair scattering and weak intra-pocket one) the gap function has 
the same sign (s-wave symmetry) within each pocket, and the nodal lines can 
happily run in between the pockets. The estimated T c for two-dimensional 
(2D) Hubbard model on such lattices is indeed almost an order of magnitude 
higher, T c ~ O.lt, as displayed in FigOHalong with the lattice structure. 

As for the dimensionality of the system, we have shown above that 2D 
systems are generally more favourable than 3D systems as far as the spin- 
fluctuation-mediated superconductivity in ordinary lattices (square, triangu- 
lar, fee, bec, etc) are concerned. Now, if one puts the idea for the disconnected 
Fermi surface on the above observation on the dimensionality, a natural ques- 
tion is: can we conceive 3D lattices having disconnected Fermi surfaces that 
have high T c 's. More specifically, can the disconnected Fermi surface over- 
come the disadvantage of 3D? If we express our idea more explicitly, what we 
have in mind is the interband nesting (or Suhl-Kondo process in its broader 
context) in the 3D disconnected Fermi surface as depicted in FigfTUI 
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Fig. 10. Interband nesting (arrow in the right panel) in 3D on a disconnected 
Fermi surface, which is exemplified here for the stacked bond-alternating lattice 
(left). 
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There, the nesting vector runs across the two bands, and this is envisaged 
to give the attractive pair-scattering interaction. So the gap function should 
be nodeless within each band, while the gap has opposite signs between the 
two bands. 

In our most recent study 109 we have found that a stacked bond-alternating 
lattice (Fig llUf l has a compact and disconnected (i.e., a pair of ball-like) Fermi 
pockets. We have shown that T c is O(0. Oli), which is the same order of that 
for the square lattice, and remarkably high for a 3D system. We have further 
found that the T c can be made even higher (~ O(0. It)) in a model in which 
the original Kuroki-Arita 2D system having disconnected Fermi surface is 
stacked. So the final message obtained here, starting from ID and ending up 
with 3D, is that 3D material with considerably high T c can be expected if we 
consider appropriate lattice structures. 

7 Closing remarks 

So we have seen the electronic properties of electron systems with short- 
range repulsive interactions, starting from ID up to 3D systems. While the 
quasi-lD ladders already contain seeds for the d-wave pairing, anisotropic 
pairing has more degrees of freedom in 2D and 3D where the topology (e.g., 
disconnected Fermi surfaces) of the Fermi surface can greatly favour higher 
Tc. Finally it would be needless to stress that the electron correlation is such a 
fascinating subject that there are many open questions to be explored. Among 
them, a question we can ask is what would happen to the superconductivity 
from the repulsive electron interaction when disorder is introduced in the 
system. Then we have a problem of dirty superconductors, i.e., an interplay 
of interaction and disorder. For ladders there are some discussions on this. 
For instance, Kimura et al. |110j have looked at the dirty double wire (i.e., 
two-band Tomonaga-Luttinger system with impurities), and noted that in the 
phases where the pairing correlation is dominant, the Anderson localisation 
is absent despite the system being quasi-lD. How this would be extended to 
higher dimensions is an interesting issue. 
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